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Abstract 

For the dilepton production at hadron collider in TeV scale gravity models, in- 
clusion of QCD corrections to NLO stabilises the cross section with respect to scale 
variations. The K-factors for the various distributions for the ADD and RS model at 
both LHC and Tevatron are presented. 

1 Introduction 

Extra dimension scenarios are now essential part of the studies of physics beyond the 
Standard Model (SM). They provide an alternate view of the hierarchy between the 
electroweak and the Planck scale. Some of these extra dimension models invoke the 
brane world scenarios to hide the extra spacial dimensions from current observation. 
Two such models that are phenomenologically widely studied are the Arkani-Hamed, 
Dimopoulos and Dvali (ADD) [1] and the Randall-Sundrum (RS) [2] models. 

In the ADD case the compactified extra dimensions could be large and the large 
volume of the compactified extra spacial dimension would account for the dilution 
of gravity in 4-dimensions and hence the hierarchy. In this model, new physics can 
appear at a mass scale of the order of a TeV. A viable mechanism to hide the extra 
spacial dimension, is to introduce a 3-brane with negligible tension and localise the 
SM particles on it. Only the graviton is allowed to propagate the full 4 + d dimen- 
sional space time. As a consequence of these assumptions, it follows from Gauss Law 
that the effective Planck scale Mp in 4-dimension is related to the 4 + d dimensional 
fundamental scale Ms through the volume of the compactified extra dimensions [1]. 
The extra dimensions are compactified on a torus of common circumference R. The 
space time is factorisable and the 4-dimensional spectrum consists of the SM confined 
to 4-dimensions and a tower of Kaluza-Klien (KK) modes of the graviton propagating 
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the full 4 + d dimensional space time. The interaction of the KK modes h^^^ with the 
SM fields localised on the 3-brane is given by 

-1 oo 

^^^nt--^Y.T'''{x)hfX^) , (1) 
^ n=0 

where T^'^ is the energy-momentum tensor of the SM fields on the 3-brane. The zero 
mode corresponds to the usual 4-dimensional massless graviton. The KK modes are 
all Mp suppressed but the high multiphcity could lead to observable effects at present 
and future colliders. 

In the RS model there is only one extra spacial dimension and the extra dimension 
is compactified to a circle of circumference 2L and further orbifolded by identifying 
points related by |/ — > —y. Two branes are placed at orbifold fixed points, |/ = with 
positive tension called the Planck brane and a second brane ai y — L with negative 
tension called the TeV brane. For a special choice of parameters, it turns out that the 
5-dimensional Einstein equations have a warped solution for < y < L with metric 
gfj,u{xf, y) = exp{—2ky) ri^i,, g^^y = and Qyy = 1. This space is not factorisable and has 
a constant negative curvature — AdS^ space-time, k is the curvature of the AdS^ space- 
time and rjni, is the usual 4-dimensional fiat Minkowski metric. In this model the mass 
scales vary with y according to the exponential warp factor. If gravity originates on the 
brane ai y = 0, TeV scales can be generated on the brane at y = L for kL ~ 10. The 
apparent hierarchy is generated by the exponential warp factor and no additional large 
hierarchies appear. The size of the extra dimension is of the order of Mp^ . Further 
it has been showed that [3] the value of kL can be stabilised without fine tuning by 
minimising the potential for the modulus field which describes the relative motion of 
the two branes. In the RS model graviton and the modulus field can propagate the 
full 5-dimensional space time while the SM is confined to the TeV brane. The 4- 
dimensional spectrum contains the KK modes, the zero mode is Mp suppressed while 
the excited modes are massive and are only TeV suppressed. The mass gap of the KK 
modes is determined by the difference of the successive zeros of the Bessel function 
Ji{x) and the scale mo = k e~'^^^ . As in the ADD case the phenomenology of the RS 
model concerns the effect of massive KK modes of the graviton, though the spectrum 
of the KK mode is quite different. In the RS model the massive KK modes h^^J{x) 
interacts with the SM fields 

Ant - -j^T^''{x)h^Kx) --rY. T''''{x)htKx) ■ (2) 

Mp Mp 

The masses of h^f^J{x) are given by Af„ = Xn k e'^'^'", where x„ are the zeros of the 
Bessel function Ji{x). In the RS model there are two parameters which are Cq = k/Mp, 
the effective coupling and Mi the mass of the first KK mode. Expect for an overall 
warp factor the Feynman rules of RS are the same as those of the ADD model. 
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These models are being tested at the Tevatron [7, 8] and will be tested at the LHC in 
the near future. One of the most important channels that could probe the signatures 
of these extra-dimensional models is large invariant mass di-lepton pair production, 
(also called Drell-Yan(DY) production) in hadronic colliders. In this process, the 
gravitons which are KK modes of the theory can appear in the intermediate state 
to produce di-lepton pairs because gravity couples to anything and everything. In 
ADD model, even though the coupling of gravity to standard model particles is small, 
the collective effect coming from all the KK modes can result in observable effect in 
some kinematic domains that can be probed. In RS, the warp factor compensates the 
small couphng making the effect feasible. Since the process under study is initiated by 
quarks and gluons, there is a large theoretical uncertainty coming from factorisation 
and renormalisation scales and also from missing higher order corrections resulting from 
strong interaction dynamics namely Quantum Chromodynamics(QCD). It is important 
for such searches to have some quantitative estimate of the effects of higher order QCD 
corrections and here we review the recent results on the NLO QCD corrections for 
various distributions of the di-lepton pair production both at the Tevatron and LHC. 

2 Di-lepton pair production at hadron collider 

Consider the collision of hadrons Pi,P2 to leptonic final states Pi{pi) + P2{P2) 
^~^{h) +^~{h) +X{Px), where X denotes the final inclusive hadronic state and carries 
the momentum Px- In the QCD improved parton model, the hadronic cross section 
can be expressed in terms of partonic cross sections convoluted with appropriate parton 
distribution functions 



The partonic level process is mediated by the SM 7 or Z boson and the graviton. 
Interestingly, the interference between the SM and graviton diagrams identically vanish 
when the phase space integration is performed [4]. To NLO the following subprocesses 
involving graviton contribute 






q + q^G + g, 



q + q ^ G + one loop , 



q + g ^G + q, 



q + g ^G + q, 



(4) 



g + g ^G + g, 



g + g ^ G + one loop . 
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The cross sections beyond LO involve the computation of one loop virtual gluon cor- 
rections and real gluon bremsstrahlung contributions to LO processes. We also include 
processes with a gluon in the initial state. Since we are dealing with energy momentum 
tensor coupled to gravity expressed in terms of renormalised fields and masses, there is 
no overall ultraviolet renormalisation required. In other words, the operator renormal- 
isation constant for the energy momentum operator is identical to unity to all orders in 
perturbation theory. But we encounter soft and coUinear divergences in our computa- 
tion. These divergences are regulated using the dimensional regularisation. We define 
n = 4: + s where n is the space-time dimension. With this procedure all divergences 
appear as l/e" where a — 1,2. The soft divergences coming from virtual gluons and 
bremsstrahlung contributions cancel exactly according to the Bloch-Nordsieck theo- 
rem. The remaining coUinear divergences are removed by mass factorisation which is 
performed in the MS scheme. The Drell-Yan coefficient function after mass factorisa- 
tion, denoted by A^^f^^z , Q"^ , /ip) , is computed by 



Ai,{z, Q\ 11 e) = rca(^, /^F, ® r,6(z, 11 e) ® A^(z, Q\ ill) , (5) 



where A^f,(^, Q^, is the bare partonic coefficient function before mass factorisation 
is carried out. The factorisation scale is given by /x^ and the kernel T^^iz) in the MS 
scheme is given by 



where = agd^ji) / 4:7r and {^) are the leading order Altarelli-Parisi splitting func- 
tions. Expanding Eq. (5) up to order a^, one can compute the coefficient function 
Al^f,{z, Q^, hf) from the bare A^^(2;, Q^, /ip, 1/e) and the known Altarelli-Parisi kernels 
P'^ ■ Finally we have to fold these finite A\^{z, Q'^, hf) with the appropriate partonic 
distribution functions to arrive at the distribution for the DY pair. For completeness 
we present the results for the beyond SM effects coming due to the graviton exchange 
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(Hqg{xi, X2, H%) + Hgq{xi,X2, AJ^^'^ (^, //|) 
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i/,,(xl,x„A.^)(A£)'«(z,Q^/.|) + a,A«■«(z,g^/.|)) ,(7) 
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the factors Tg correspond to pure gravity part, Rah are the renormahsed parton dis- 
tributions and l^ah^ are given in [4]. 



3 Discussion 

The NLO-QCD corrections have been calculated in the ADD case for various distribu- 
tions and for forward-backward asymmetry of the invariant lepton pair at both LHC 
and Tevatron [4]. It was further extended to the RS model in [5]. The un-integrated 
distribution with respect to cos Q* (the angle of a lepton with respect to the hadron in 
the cm frame of the lepton pair) is particularly important as the interference between 
the SM and gravity is not zero any more [6]. The inclusion of the higher order QCD 




Figure 1: The various contributing sub-process to the invariant mass distribution of 
the dilepton pair for ADD model at LHC and Tevatron. 

corrections reduces the dependence on the factorisation scale significantly making the 
theoretical predictions for the relevant observable more reliable for physics studies at 
high energies in hadron colliders. The dominant production mechanism for di-lepton 
production is through gluons in the initial state in the case of LHC and quarks in 
the case of Tevatron. The parton distribution functions corresponding to these initial 
states at LHC energies usually have large factorisation scale dependence. We find that 
our results on NLO QCD corrections to the observable mentioned in this paper do re- 
duce the factorisation scale dependence coming from the parton distribution functions 
as expected. We also find that the NLO QCD corrections quantified by the K factor 
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are not small. The K-factor at large invariant mass is found to be as large as 1.6 at 
the LHC while at the Tevatron, it is similar to the SM value [4]. This large value at 
LHC is again due to the initial state gluons which have large flux at LHC energies. 
The behavior of the K-factor for RS model can be understood from the fact that only 
in the resonance region the gravity part contributes, so at off resonance it is the SM 
K-factor while on resonance it is dominated by the gravity part. 
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